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Abstract 

We present new evidence in support of the Penrose's strong cosmic 
censorship conjecture in the class of Gowdy spacetimes with T 3 spa- 
tial topology. Solving Einstein's equations perturbatively to all orders 
we show that asymptotically close to the boundary of the maximal 
Cauchy development the dominant term in the expansion gives rise to 
curvature singularity for almost all initial data. The dominant term, 
which we call the "geodesic loop solution" , is a solution of the Ein- 
stein's equations with all space derivatives dropped. We also describe 
the extent to which our perturbative results can be rigorously justified. 
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1 Introduction 



It is well known p]] that in the class of globally hyperbolic spacetimes arbi- 
trary Cauchy data for the vacuum Einstein equations have unique maximal 
developments (maximal Cauchy developments). We would like to know un- 
der which conditions a similar statement could be true in some larger class of 
spacetimes, possibly all vacuum Einstein spacetimes satisfying certain differ- 
entiability conditions (see []|] for a nice review). The Taub-NUT spacetime 
H is prototype of a spacetime where the globally hyperbolic part (Taub) can 
be extended into at least two non-isometric Taub-NUT spacetimes ffl, thus 
violating the uniqueness of the maximal development. Penrose's Strong Cos- 
mic Censorship (SCC) conjecture states that for almost all (generic) initial 
data the maximal development is, indeed, globally hyperbolic, not larger, 
hence restoring the uniqueness of the maximal development. The spacetimes 
like Taub-NUT obviously necessitate the exclusion of some spacetimes as 
special (nongeneric) . In the case of Taub-NUT it is the high symmetry (its 
isometry group is four dimensional) that makes it nongeneric. 

To prove the SCC we would have to specify arbitrary Cauchy data, find 
the maximal globally hyprebolic development, and show that for almost all 
initial data that development is inextendible by proving, for example, that: 

1. The maximal globally hyprebolic development is geodesically complete 
or 

2. For every incomplete geodesies some curvature scalar blows up when 
approaching the incomplete end or ends (the singularity). 

The Hawking- Penrose singularity theorems tell us that, generically, the 
first case is not true, but they do not tell us anything about the behav- 
ior of the curvature when we approach the incomplete end of a geodesic 
(the singularity). This is not surprising since these theorems use Einstein's 
equations to a very limited extent. To find the behavior of the curvature 
when approaching the singularity we'll have to find the detailed asymptotic 
behavior of the solutions of Einstein's equations close to the singularity. 

The proof of the SCC for the general case is presently out of our reach 
due to the complexity of Einstein's equations, but, to gain some insight into 
the general case, we can try to prove the SCC in the restricted class of 
spacetimes having an n-dimensional isometry group, with least possible n. 
Following that strategy in this paper we analyse the asymptotic behavior 
of the solutions of Einstein's equations in the class of Gowdy spacetimes on 
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T 3 x R; spacetimes that have 2-dimensional U(l) x U(l) isometry group 

We develop a new method for studying the asymptotic, singular behavior 
of large classes of cosmological solutions of Einstein's equations. The method 
consists principally in showing that one can solve the n-th order perturbed 
Einstein's equations for all n using the "variation of constants" method, 
and then use the freedom in choice of the n-th order solution to obtain an 
asymptotic (when approaching the singularity) sequence of functions which 
is presumably the asymptotic expansion of some exact solution of Einstein's 
equations. This sequence is uniquely determined by the zeroth order solution 
which we call the "geodesic loop solution". The zeroth order equations are 
just the Einstein's field equations with "space" derivatives dropped and the 
method is applicable only to the spacetimes whose dynamics is "velocity- 
dominated" in the sense of Eardley, Liang and Sachs (8|, || . 

For the first two orders our method gives the same results as the method 
suggested by Cosgrove, based on the "multiple scales" technique of applied 
analysis [10]. The full extent of applicability of our method is not yet known. 
In the present paper we develop and apply it to the case of Gowdy metrics on 
T 3 x R whose asymptotic behavior has been analysed before by Mansfield, 
using a different perturbative method jn|]. In a separate work we shall 
apply it to U (1) symmetric spacetimes. We believe, but this has not yet 
been explicitly demonstrated, that our basic method should be applicable 
to general, non-symmetric solutions of Einstein's equations. 

An advantage to starting with with the Gowdy metrics on T 3 x R is 
that much is known about their behavior. In particular on can prove a 
"global existence theorem" to the effect that, in a suitably rigid coordinate 
system (in which the time coordinate r measures the geometrical area of the 
two-tori which arise as orbits of the isometry group), all sufficiently smooth 
Gowdy solutions extend globally to all r £ (0, oo) and that this interval of 
existence exhausts the maximal globally hyperbolic development ||12|| . The 
same reference shows that the r = const, hypersurfaces always approach 
a "crushing singularity" of uniformly diverging mean curvature as r — > + 
but approach a boundary of infinite three volume as r — > oo. The former 
boundary always arises after a finite lapse of proper time ( as measured, 
say, along the normal trajectories to the chosen foliation from some non- 
singular reference Cauchy surface) making the spacetime past geodesically 
incomplete, whereas the latter occurs only after an infinite lapse of proper 
time. Our main interest here is with the "crushing singular" boundary which 
occurs as r — > + , we would like to show that the spacetime is inextendible 
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beyond that boundary . It is expected to be curvature singular generically, 
but may in some special cases instead correspond to Cauchy horizon across 



which (analytic) extensions are possible |T3|, 14, |15| . 

The mathematical methods used in proving the main results of reference 
[12 1 give only rather weak information about the behaviors of the solutions 
near their crushing boundaries at r = 0. For the special case of polarized 
Gowdy metrics (which are governed by a linear hyperbolic equation) on can 
strengthen the arguments of Ref . |T^] to derive genuinely sharp estimates for 
the behaviors of the metric functions and their derivatives near the crushing 
boundaries ||. In this case on can rigorously characterize the asymptotic 
behavior of every solution and even classify the solutions in a natural way by 
studying the asymptotic behaviors of their curvature tensors. While we have 
considerable hope that similar arguments can be developed for the general, 
non-linear Gowdy metrics, this has not yet been accomplished. 

Lacking a direct (i.e., non-perturbative) means for establishing the asymp- 
totic, singular behavior of the generic Gowdy spacetime, we have developed 
the perturbative approach to be presented here. This perturbative approach 
has (as we have mentioned) the advantage of applicability far beyond the 
rather limited class of Gowdy metrics. When restricted to the Gowdy class 
however, it yields results which may eventually be provable by direct, non- 
perturbative methods. For the larger classes of spacetimes to which the per- 
turbative method definitely applies (e.g., those having one spacelike Killing 
field), the possibility of a direct (non-perturbative) proof of the correspond- 
ing results seems currently to be rather remote. Thus the Gowdy metrics 
provide a natural test-case for the perturbative method. The perturbative 
results provide natural conjectures which one can reasonably hope to prove 
by direct methods and the direct methods, if successful, can be expected to 
justify the (less rigorously founded) perturbative approach. 

The outline of the rest of the paper is as follows: In section 2 we define 
the Gowdy spacetimes and briefly review the field equations. In section 3.1 
we formulate the perturbative method, in section 3.2 we solve the zeroth and 
first order equations and find the asymptotic behavior of the corresponding 
solutions. In section 3.3 we show that, as r — > + , all higher order solutions 
decay faster than the zeroth order solution, called the "geodesic loop solu- 
tion". In section 3.4 we show that for almost all initial data the "geodesic 
loop spacetimes" , which presumably capture the asymptotic behaviors of the 
exact solutions of the Gowdy field equations, are curvature-singular when, 
r — > + , and hence inextendible beyond their maximal Cauchy develop- 
ments. In section 4 we describe known rigorous results, which all support 
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our perturbative analysis. 



2 Gowdy T 3 x R Spacetimes 

Consider a spacetime that can be foliated by a family of compact, connected 
and orientable spacelike hypersurfaces. If the maximal isometry group of 
the spacetime is two-dimensional and if it acts invariantly and effectively on 
the foliation than the isometry group must be U(l) x U(l). Moreover, the 
foliation surfaces must be homeomorphic to T 3 , S 1 x S 2 , S s or to a manifold 
covered by one of these, and the action of the U(l) x U(l) is unique up to 
a diffeomorphism. In such a spacetime the Killing vector fields K a , a = 1, 2 
associated with the isometry group must commute, and two scalar functions 

c a = e af3y5 KfK^K s a (1) 

must be constant [16]. The spacetimes which satisfy the above mentioned 



symmetry requirements, and in which both constants c a vanish, will be called 
Gowdy spacetimes 0. 

For the simplest case of T 3 spatial topology, Gowdy ]fjj showed that in 
such a spacetime a coordinate system can be found such that the metric can 
be expressed in the form 

ds 2 = e 2A [-e~ 2t dt 2 + d6 2 ] + e~* [cosh W + cos $ sinh TV] (fix 1 ) 2 (2) 
+ 2e~*sin$ sinh W dxdy + e~* [cosh W - cos <J> sinh W](dx 2 ) 2 , 

where d/dx 1 and d/dx 2 are commuting Killing vector fields on T 3 in the 
S 1 x S 1 directions, and the functions W , $ and A depend only on the 
two remaining coordinates: "time" t and "angle" 9 parametrizing the third 
S 1 factor in T 3 . For convenience, instead of the time r that measures the 
geometric area of the two-dimensional group orbits we use the time t = 
— lnr. The singularity at r = 0, mentioned in the introduction, corresponds 
to t — > oo. 

Einstein's vacuum field equations give the equations of motion for the 
metric functions W, $ and A 

W, u ~ \ sinh W*J = e~ 2t [W> - i sinh 2W <S> 2 ] (3) 

$ it + 2 coth W Wj$,t = e' 2t + 2 coth W Wfi$,o] (4) 

A )t = -[1 - W 2 t - sinh 2 W $ 2 t - e- 2t {W 2 e + sinh 2 W <5> 2 e )} (5) 

= -\[W, t W fi + $, t $, e sinh 2 (6) 
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which have to be solved on the t9 cylinder. 

To solve this system of equations we proceed in two steps. We first solve 
equations (||) and (Q) for W and <£. These two equations can be put into 
Hamiltonian form with the Hamiltonian 

H[X,n] = / Hd9 (7) 





H = l -[g AB Ii A Ii B + e~ 2t g AB X A 'x B 'l (8) 

where X 1 and X 2 are W and <3? respectively; prime denotes differentiation 
with respect to 6 and g AB is the metric of two-dimensional hyperbolic space 

g = dW 2 + smh 2 Wd<S> 2 . (9) 

From Hamilton's equations it follows that imposing periodicity in 9 (with 
period 2ir) on initial conditions is enough to obtain W and $ that are 
periodic at all times. A solution W(t,9) and &(t,9) can be conveniently 
pictured as a closed loop moving in the two-dimensional hyperbolic space 
with speed 



/g AB U A U B . (10) 
After obtaining W and <3? we can find A by evaluating the line integral 

A(t,9) = a + J^a (11) 

over some path T on the tO cylinder with the end points at (to, 6q) and (t, 9), 
and with components at and ag of the two dimensional one-form a equal to 
the right-hand sides of the equations (||) and (||) respectively. Expressed in 
Hamiltonian variables 

ch = (12) 

a e = -\x A 'll A . (13) 

As a consequence of the form of the Hamiltonian the function ag satisfies 
the identity 

da e fSH\' dag 8H 1 , 

ax* IotiJ " drT A Jx* ~ ~2 n ~ at > e (14) 
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for arbitrary functions X A and n^. For X A and IT^ that satisfy the Hamil- 
ton's equations of motion the left-hand side of ( |i~4|) becomes the time deriva- 
tive of ceg and we obtain 

ote,t = OL t ,e or da = 0, (15) 

which guarantees that the line integral does not depend on the path. To 
have the function A globally defined on the t6 cylinder we must impose one 
more constraint on W and namely that A must also be periodic in 9 with 
period 2tt, which leads to the constraint 

a e {t,Q)dO = Q. (16) 

o 

This constraint has to be imposed only at one time to because it is conserved 



during the time evolution due to equation (|15|) . 

Two important subclasses of Gowdy spacetimes are obtained by impos- 
ing special conditions that are conserved by the evolution equations. We 
shall refer to the spacetimes satisfying the condition = as the "polarized" 
Gowdy spacetimes and to the spacetimes satisfying W = W(t) and $ = n9, 
for some n € Z, as the "circular loop" spacetimes. For the polarized case 
the evolution equations for W and $ reduce to one linear partial differential 
equation, and for the circular loop case to one non-linear ordinary differential 
equation. 

To facilitate further analysis we will introduce new coordinates x and y 
in the hyperbolic space, such that 



tanhw = VSrtw (17) 

tan $ = x2 + j/2 ~ 1 , (18) 

These are just the Poincare half-plane coordinates in which the hyperbolic 
metric becomes 

dx 2 + dy 2 „ . . 

g = x£R, y>0. 19 

The equations of motion for x and y are also given by Hamilton's equations 
where now, however, the metric in (|8|) is replaced by expression ( |l9| ) and 
{X 1 , X 2 } now stand for {x, y}. With this modification the Hamiltonian den- 
sity becomes a rational function of its arguments. The functions appearing 
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in the Gowdy metric (||), expressed in terms of the x and y, are 
cosh W + cos sinh W 
cosh W — cos sinh W 
sin <3? sinh W 

Note that the y > restriction is necessary to obtain a metric (0) that 
is of Lorentzian signature and that the transformation x — > — x gives an 
isometric metric ( it amounts to exchange of x 1 and x 2 ). 

3 The Perturbative Expansion 
3.1 The Perturbative Method 

As already mentioned, in general we do not know how to solve the evolu- 
tion equations for x and y. What we can do, however, is try to solve the 
equations perturbatively to all orders and analyse the asymptotic behavior 
of the expansion to learn about the behavior of the Gowdy space-times close 
to the singularity at t = oo. To facilitate this analysis we will first intro- 
duce a seemingly even more complicated problem. We will try to solve the 
evolution equations for the e dependent Hamiltonian density 

H e = \[g AB U A U B + e e- 2t g AB X A 'x B '}. (21) 

The case of real interest to us is e = 1, but it is easy to see that any 
solution of the e-dependent equations for some ei > is also a solution 
of the equations for any other e 2 > , provided we shift the time t by an 
appropriate constant. Precisely stated, if X ei (t) is a solution of the equations 
of motion for e = e\ > then 

X e2 (t)=X tl (t-hne 2 /e 1 ) (22) 

is a solution of the equations for e = e 2 > 0. So, with the introduction of e 
we have not changed any significant property of the problem, but the e in the 
Hamiltonian will help us obtain simple recursive form for the perturbative 
calculations. 



(x + y) 2 + 1 
{x-y) 2 + l 



2y 



x 2 + y 2 



1 



2y 



(20) 
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Assume that there exists an e-dependent family X(e,t,9) of solutions 
to the modified equations which depends smoothly upon the parameter , 
in particular around e = 0, and expand it into a power series in e to find 
evolution equations for each term in the expansion: 



X(e,t,6) 
X (n \t,d) 



e=0 



X 



n=0 

t-X(M,0) 

(*\ 

y 

X 

V V J 



(23) 
(24) 

(25) 



Whether the family actually exists depends on the convergence properties 
of the "formal" perturbative expansion ( ]23| ) which can be changed using 
the freedom in choosing appropriate X^ in every order since the n-th or- 
der perturbation equations don't have unique solutions; but as we shall 
see, irrespective of the convergence of the expansion, the recursive calcula- 
tions of the higher order terms are always well defined in terms of the lower 
order terms. It could easily happen that the expansion could not be made 
convergent but could nevertheless be made asymptotic as t — > oo, in which 
case we would still be able to study the behavior of Gowdy spacetime close 
to the singularity. The question still remains whether we can study all, or 
almost all, Gowdy spacetimes using this perturbative method; i.e. can we 
for any Gowdy solution find an e-dependent family which for e = 1 reduces 
to that particular Gowdy solution and is smooth around e = 0? We will 
return to this question later, for now let us just say that a function counting 
argument, together with the properties of our expansion, suggests that the 
expansion is possible for an open subset of too strong all Gowdy solutions 
on T 3 x R. 

To clarify the analysis let us for a moment consider an s-dimensional 
system of nonlinear partial differential equations on R x M, with "time" 
coordinate in R and "space" coordinates in an arbitrary finite dimensional 
manifold M. In shorthand notation we write 



X = h(e,X,Y)=g(X) + ef(Y), 



(26) 



where X is an s-component function and Y is a multicomponent function 
that contains all the components of X and some number of partial derivatives 
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of X with respect to the "spatial" coordinates. For e = this system reduces 
to a system of ordinary differential equations. By differentiating ( |26|) with 
respect to e and then putting e = we obtain the equations that have to be 
satisfied by if a family X(e) exists: 



d" 



*W = ^Ke,X,Y) 



h (n) 



6=0 



n = 
n > 



xM=g(xW) 

X(n) = g (n) + n f(n-l) 



(27) 
(28) 

(29) 



Now we will analyse the above perturbative equations without worrying 
about the existence of an e-dependent family of solutions of the equations 
(|2~6|). The derivatives of g(X) and f(Y) can be calculated using the chain 
rule, and with all the indices explicitly written, the expression for g{X) is 
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(n) 



e=0 



fc=l ii— ijt 



(mi) 



a fc ^ r ., fc (x(°)) 



• • • dXi 



(30) 



(31) 



e=0 



From now on we will suppress the vector indices, like ji\ • • - ifc, and argu- 
ments of d n g and d n f to simplify the notation. For the same reason we 
will put g(X^) = #(°) and f(Y^) = /(°). In this shorthand notation the 
evolution equations for X^ for any n > 1 are 



fc=2 



m\---m k 



n-1 



k=l 



y(mi) . . . y(m fe ) 



^m fc =n-l 



(32) 



The first three equations including the zeroth are 
n = X(°)= 5 (XW) 



(33) 
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n = l iW= dgX^ + f(Y^) (34) 
n = 2 = dgX® + d 2 gX^X^ + 2dfY^ (35) 

We see that to all orders, except the zeroth, the perturbative equations 
are ordinary inhomogeneous linear differential equations of the form 

X {n ^ = dg X^ + S {n \ (36) 

where the matrix that multiplies X^ is the same for all orders and the 
source depends only on the lower order functions X^ and 7", i.e. 
k,l < n. Once we know the general solution for the zeroth order equation 
(|33|), which we shall refer to as the seed solution, we can recursively solve 
all the higher order equations using the "variation of constants" method for 
ordinary linear differential equations. The general solution of the n-th order 
equation (^) is 

X {n) (t) = Xj^jt) + fg{t,t') S^(t')dt', (37) 



where xj™l, is a solution of the homogeneous part of Eq.(36) and Q is its 



horn 

Green's function. Both can be calculated from the seed solution. By differ- 
entiating the seed solution with respect to the constants c, it depends on, 
and linearly combining these derivatives, we obtain the general homogeneous 
solution 

da 



■V 



Xhom = X! ai ~ ~„ — ' (^) 



where on are arbitrary constants. The Green's function Q(t,t') for (|37|) is 
the Jacobian matrix for the transformation 

gggff (39) 

and can be calculated from the seed solution expressed in terms of the initial 
values. It is easy to check that 

Xhom = 9g X hom (40) 

and 

^G{t,t') = dg{t)Q{t,t') W, (41) 
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and confirm that (^7() satisfies (pq). The lower limit of integration in < pT\ ) 
could have been chosen to be different at every order as well, but a change in 
its value corresponds just to adding a solution of the homogeneous equation. 

Using any finite number of the solutions of the perturbative equa- 
tions (^7|) we can construct an e-dependent function 

N n 

which we shall refer to as the formal expansion of X . For any formal 
expansion we can define a function 

5 N (e)=X N (e)-h(e,X N ), (43) 

which, as a consequence of perturbative equations, satisfies 

^Me)U^4 n) = 0, Vn < N. (44) 

To simplify the notation we have assumed that h depends only on X. Given 
a formal expansion X/v of X, we can generate a formal expansion Fn of any 
smooth function F by 



= E - { F{n) i F(n) = ^nx N {e)) 



n=0 



(45) 

e=0 



where F^ in general depends on all X^ k \ with k < n, and is explicitly 
given by (^) with g replaced by F. 

An interesting consequence of the above definition and Eq.(|43) is that 
if we have a constant of motion for the evolution governed by (|2q), i.e. a 
function C{X) such that identically for any X and e 

^h(e,X)=0, (46) 

all the terms C'"' in the formal expansion of C will be constant in time 
also. To see that, consider a formal expansion CV of C. Using (46) we can 
evaluate 

C(X N ) = d -^X N = ^l [ h(e,X N ) + S N (e)] 

dC(X N ) ' 
= q x Me), (47) 
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differentiating n times with respect to e and using (44) we obtain 



= 0. (48) 

e=0 



The general idea of using 5 N = could be employed to find other properties 
of the perturbative equations that are inherited from the "full" equations. 

To conclude this general discussion we once more note that one has con- 
siderable freedom in choosing the formal expansion apart from the freedom 
in the choice of the seed solution. One can use this freedom to try to obtain 
an expansion in which, when time goes to some to, all terms in the expansion 
form an asymptotic sequence of functions i.e. 

. w \ ; ^ t -» t , (49) 

for all n; or perhaps even to obtain a convergent expansion. 

Since we are interested in finding the asymptotic behavior of the solutions 
of the Gowdy equations when t — > oo, we shall chose the lower limit of 
integration in (]37|) to be to = co and let X^ m = for all n. This choice leads 
uniquely to the fastest possible decay of the higher order terms, because the 
homogeneous solutions of the perturbative Gowdy equations do not decay 
faster than the seed solution when t — > oo. The in (^) after shifting 
the integration variable t' — > t' — t becomes 

POO 

(t) = - g(t,t + 1') sw (t + t')dt', (so) 

Jo 

and depends only on the seed solution. As we shall see, this choice leads to 
a uniformly behaved asymptotic sequence for a large range of seed solutions, 
which suggests that the sequence might be an asymptotic expansion of an 
exact Gowdy solution. 



12 



3.2 The Geodesic Loop Solution and the First Order Cor- 
rection 



Returning to the Gowdy equations we note that Hamilton's equations for 
the Hamiltonian (^) are of the form ( 26 ) with 



X 





t 


X 


\ 


( x \ 






y 




y 


Y = 




x' 




X 






y' 




V y J 






x" 






\ 


y" 


) 



(51) 



and 



9(X) 



where 



^3 

x 4 

A 

\C-BJ 



( 



fOO 



-2t 



o 
o 



D-E 
\F+G-H j 



A 
B 



X 2 



X3X4 



D 



-Xi 



C =x^ 



E = —Y 3 Y 4 
Y 2 

F = Y 6 

G = —Yi 
Y l ^ 

H = —Y? 

Y 2 4 



(52) 



(53) 



The zeroth order equations are, as is obvious from the Hamiltonian (21), 
just the geodesic equations in the Poincare half-plane whose general solu- 
tion, using our notation, is 



X 



(0) 



X. 



(0) 



a + &tanh(ct + d) 
1 



X. 



(o) 



be 



cosh(ct + d) 
1 



X. 



(0) 



cosh 2 (ct + d) 
sinh(ct + d) 



-be 



cosh (ct + d) 



(54) 
(55) 

(56) 
(57) 
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where a, b > 0, c > and d are time independent but otherwise arbi- 
trary smooth functions of 9. For convenience we shall refer to these as 
the "geodesic loop solutions", since each curve I™(t, 9q), for fixed 6>o, is a 
geodesic. The speed of these geodesic loops, evaluated using Eq. (|i~0|) , is 
equal to the function c(9) and is, as it should be, constant in time. The 
geodesic loop solutions have the same number of free initial data as the gen- 
eral Gowdy solutions. Note that when a, b, c and d do not depend on 9, the 
geodesic loop, in this case degenerated to just one point, is a solution to the 
full Gowdy equations. The geodesic loop solutions give rise to the "geodesic 
loop spacetimes" which, we hope to prove, are asymptotically (as t — > oo) 
approached by Gowdy spacetimes. 

Differentiating the geodesic loop solutions with respect to a, b, c and d 
we obtain the general homogeneous solution for the perturbative Gowdy 
equations. It is easy to see that it does not decay faster than the seed 
solution when t — > oo, which justifies our choice of in (|50| ) as unique 
for a (potentially) well behaved asymptotic expansion. 

Since we are interested in the asymptotic behavior when t — > oo we will 
expand the above functions in power series in l/£, where £ = e^ ct+d \ which 
are convergent for all £ > 1 i.e. for all t and 9 for which ct + d > 0. To 
exhibit the asymptotic behavior we will factor out the lowest power of £. 
For example: 

oo 

X[ 0) = a + 6+^26(-l) fc r 2fc (58) 

k=l 

oo 

xjp = r^K-r/r 2 *- (59) 

k=0 

In fact all the the zeroth order functions have a similar form. The other two 
components of X^ are time derivatives of the first two, and includes 
also 9 derivatives of X^ . 

Definition 1 Let F:R — > R be a function that can be expanded into a uni- 
formly convergent series 

oo 

F(t)=C N Y,akC* k , (60) 

k=0 

where £ = , N G R , all ak are polynomials in t of degree less than or 

equal to some fixed integer p > and oq ^ 0. Then, for easy reference, we 
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call F a ^-expandable function , N(F) the decay exponent of F , p(F) the 
polynomial degree of F and ao(F) the dominant coefficient of F. If F(t) = 
for all t we put N(F) 



oo. 



This definition will be important in the inductive proof of the decay of all 
higher order terms since all the functions we will be dealing with in this paper 
satisfy the conditions of the above definition, i.e. they are ^-expandable. For 
any two such functions F and G the product FG is obviously a ^-expandable 
function, with p{FG) = p{F) + p{G) and with N(FG) = N(F) + N(G). 
Moreover, if N(F) — N(G) = 2k for some integer k > the sum F + G 
will also be ^-expandable and N(F + G) = N(G). In the special case 
when N(F) = N(G) we have N(F + G) > N(F), because of the possible 
cancellation of the dominant coefficients. If k dominant coefficients cancel 
we have N(F + G) = N(F) + 2k. So, in general, for any two ^-expandable 
functions F and G whose sum is also ^-expandable, we have 



N(F + G)> mm{N(F),N(G)} 
Owing to the uniform convergence of expansion (|60[) the integral 



(61) 



/ F(t)dt = / £ a k r 2k ~ N dt = £ a k r 2k ~ N dt, (62) 
Jo Jo k=Q 



k=0 



for all F for which N(F) > 0. Therefore, to evaluate ( |50| ) we need to 
integrate only simple integrals of the form 



fe -{2k+N)ct dt 



[{2k + N)c 



r+l 



(63) 



where the integer k > and the integer r < p(F), the polynomial degree of 
F. Note that all the zeroth order functions X^ ' have p(X^) = 0, and that 
each differentiation with respect to 9 brings one t down from the exponent 
and therefore increases the polynomial degree by 1. 

In the generic case, i.e. when no special conditions are imposed on the 
functions a, b, c and d, the decay exponents for the zeroth order functions 
are 

( o \ 



N(X 



(oh 



1 

2 

V i J 



N(Y 



(Oh 



1 

1 


V 1 J 



(64) 
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from which follows that 



N(A) = 2 
N(B) = 3 
N(C) = 1 



N(D) = 
N(E) = 
N(F) = 1 
N(G) = -1 
N(H) = 1. 



(65) 



From (34) and (52) we can calculate the first order source vector 
its decay exponents are 



N(S {1) ) = N(fW) 



I oo 

oo 

A + 2 
V A + l / 



(66) 



where A = 2(1 — c)/c, and its polynomial degree p(S^) = 2. 

Prom the form of Gowdy's g and / we see that the first two components of 
the 4-dimensional source vector are zero for all k, which implies that we 
need just the last two columns of the Green's function to evaluate integrals 
(|50|). Explicit calculation of the Green's function 4x4 matrix Q(t,t + t') 
shows that it depends on 8 only through the functions c and d, and that 
all its components are of the form ( |60| ) with the coefficients polynomials 
of the first order in t' multiplied by e nt , with n ranging from -1 to 2. The 
time t appears only in the exponents, therefore the polynomial degree of Q 
is zero. Asymptotically for t — > oo we obtain 



V 



(l-C 2 )/2c 

r 1 (i-c 2 + 2ct / )/c 
c 2 

rH3C 2 -2ct'-3) 



r 1 (r 1 -c + 2Ccf)/c \ 
aa' + 1) / 



(67) 



Here we have put asterisks for uninteresting components and we have written 
( for e ct . The highest power of £ in the components of the Green's function 
will be important for the convergence properties of the integral (|50|), and 
it is convenient to introduce another simple definition to keep track of the 
exponential functions of the integration variable t' . Let N'(F) denote the 
function defined like N(F) in Definition 1, with t changed to t' and with 



£ = e ct> written in place of £ = e 



(ct+d) 



All the functions we use obviously 



satisfy the conditions stated in Definition 1, with respect to both t and t'. 
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Pairing the exponents N and N' in parentheses for brevity, from the exact 
expression for the Green's function we obtain 



( * * (0, 


-2) 


(1, 


- 1} \ 


* * (1, 


-2) 


(0, 


-1) 


* * (0, 


-2) 


(1, 


-1) 


\ * * (1, 


-2) 


(0, 


-1) ) 



(68) 



Because the identity N(F(t + t')) = N'(F(t + t')) is valid for any function 
satisfying the conditions in Definition 1, using (|66[), we obtain 



(N,N')(SV>) 



( (oo,oo) 
(oo, oo) 
(A + 2, A + 2) 
V (A + 1,A + 1) J 



(69) 



Multiplying the Green's function and the source we find the integrand whose 
exponents are 

( (A + 2,A) \ 



(A + 1,A) 
(A + 2, A) 
V (A + 1,A) J 

The integral will converge for all for which 



A(0) >0^c(6) € (0,1) 



(70) 



(71) 



Finally, the integration gives the first order functions which are, because of 
(|62"D, also ^-expandable. We were able to calculate the first order integral 
explicitly, without restoring to the ^-expansion, and thus obtain the exact 
expressions for the first order functions, whose asymptotic form agreed with 
the result obtained by ^-expansion. The polynomial degree of JfW is 2, the 
same as p(S^), since p(G) = 0, but the dominant coefficient of all lW-s 
are just first degree polynomials in t. The differences of the decay exponents 
of the first and zeroth order functions are 



N(XM)-N(XW) 



(o)> 



/ A + 2 \ 
A 
A 

V A J 



N(Y 



N(Y 



(0)^ 



/ A + 2 \ 
A 

A + 2 
A 

A + 2 
V A J 



(72) 
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and, again, A > is the condition that has to be satisfied if we want that 



X 



(i) 



X 



(0) 



te 



-2(l-c)t 



when t — > oo, 



(73) 



for i = 1,2, 3, ( x[^ decays even faster). 

The condition c{9) < 1 can be relaxed for 9 in some Iq C [0,27r), if 
a{9) + b{9) =const. for 9 G Iq. Explicit calculation in this case shows that 
the decay exponents of the first order source increase by 2 and the first order 
integral converges for all c{6) > , 9 S Iq giving 



(oh 



/ A + 4 \ 
A + 2 
A + 2 



(0)^ 



/ A + 4 \ 
A + 2 
A + 4 
A + 2 
A + 4 
V A + 2 J 



(74) 

Now when we know the decay exponents of the zeroth and first order 
functions we can go on to calculate the decay exponents of the higher or- 
der functions using d5(]), but all we actually need to know is whether they 
increase with order. 



3.3 Higher Order Terms 

The lower bounds for the decay exponents of X^ can be easily found for 
any order m following a rather simple inductive argument. 

Theorem 1 Let be a solution of the perturbative Gowdy equations 

given by ^5^), and let its geodesic loop parameters a, b, c and d be arbitrary 
smooth functions of 9 except that c{9) € (0, 1), then for any m > 1, the 
decay exponent of X^ m > satisfies the inequality 

N(X^)>N(X^) + (m-l)X. (75) 

Proof: The statement of the theorem is trivially true for m = 1. 

Inductive hypothesis: Assume that the statement ( |75| ) is true for all 
1 < m < n — 1. Using the result (|72|) from the previous section we obtain 

N(X^) > N(X<®) + mX, (76) 
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for all 1 < m < n — 1. Now, all we need to prove the theorem is show that 
the decay exponents of the n-th order source satisfy the inequality 

iV(5 (n) ) >iV(5 (1) ) + (n-l)A (77) 

which will, because of the inequality (|6l|), after multiplication with the 
Green's function and integration over t' give the desired result for 
m = n. 

The components of the source vector for n > 1 are sums of parts of 
the form 

S {n) (g) + S (n) (/) (78) 

S (n) ^ = Qk g X (mi) . . . X (m h ) ^ 

s^ n \f) = d k f y (mi) • • • y(™<=) (so) 

22 m;=n-l 

where for g we can substitute A, B or C and for f D, E, F, G and H. 
Again, because of the inequality (|6l|), we need to find only the minimal 
decay exponent of any term in the sum to obtain a lower bound to the 
decay exponent of . To do this we need to take a closer look at the form 
of the terms above. 

All the functions A, B, C, D, E, F, G and H are of the form P = UjZf , 
where Z stands for X and Y, e.g. A = 2X^X3X4. Derivatives of such 
functions are easy to evaluate and we have 

f)P k- f)P 

- - - — P => N(-j£r) = N(P) - N(Xj) (81) 



dXj Xj dXj 



if kj ^ 0, otherwise N = oo. Using the above property we find 
N(d k f y (mi) • • .y( m *)) = A(/ (0) ) 



k 



+ [A(y (m °) - A(y^)] (82) 



1=1 

2~2 mj=n-l 



N{d k g X^ ---X^) = N(gW) 

k 

+ [N(X^)-N(X^)] (83) 



i=l 
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In the sums above only X-s and Y-s of the order n — 1 or lower appear and 
the inequality (^) together with the conditions J2 m i = n — 1 and J2 m i = n : 
for Y and X sums respectively, enable us to conclude that the sums satisfy 
the following inequalities 

k 

^[jV(y(«w))_jV(y(o))] > (n-l)A, (84) 

i=i 

fc 

J2[N(X^) - N(X^)] > nX. (85) 

Recalling that = /(°) we obtain 

N(S( n \f)) > N(S^) + (n-l)X, (86) 
and using N(g^) = JV(/( )) - A and fl8|) we get 

N(S^(g)) >N{sW) + (n-l)X. (87) 

Therefore, the source S^"-* has the desired behavior, which concludes the 
proof. 

When no cancellations in the dominant coefficients are present in the 
course of calculations, the statement of the theorem becomes an equality, 
and for n > 1 and for t — > oo we obtain 

x (n) = x W, &e -Mi-c)^ (88) 

yW = } ~ q 2n e- 2n ^ t e- ct . (89) 

The p2n and q2 n are polynomials in t of degree 2n or possibly lower, which 
can be easily shown using an argument similar to the one used for the decay 
exponents. We see that the terms decay faster and faster with order and that 
we have obtained an asymptotic sequence. When cancellations are present, 
the decay is even faster and we can group the terms to obtain an asymptotic 
sequence again. 

A statement similar to the Theorem 1. is true for the case a + 6=const. 
and c arbitrarily large. Then the higher order terms decay even faster, and 
their decay exponents satisfy the inequality 

N(X^)>N(X^) + (n-l)(X + 2) (90) 
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3.4 Expansion of the Metric Functions and Curvature 

Having obtained the asymptotic sequence of solutions to the perturbative 
equations for the Gowdy functions x and y we can generate the expansion 
for any smooth function of x and y using fl45|). For example, let F be a 
function of X of the form 

f « = Wi' <91) 

where P and Q are polynomial functions of X. To find the expansion 

M n 

F M = m (0) ) + Ei F(n) ( 92 ) 



, n! 

n=l 



to arbitrary finite order M, we have to calculate the F^ n '-s. They are given 
in terms of d k F by the equation (^), with g replaced by F. The first partial 
derivative of F, 

~dP_ dQ~ 
.P ~ Q . 



8F=P Q 



(93) 



is again a function of the form P/Q, with P and Q polynomials in X. For 
each monomial T in, say, P we have N(dT) = N(T) — N(X), which using 
(fip gives N(dP) > N(P)-N(X). Since N(dP/P) =N(dP) - N(P) > -N(X), 
using we obtain 

N(dF) > N(F) - N(X), (94) 

and repeating the same reasoning k times, 

N(d k F) > N(F) - kN(X). (95) 

Using the definition of F^ and Theorem 1 together with ( |95|) we obtain 
that higher order terms in the expansion of F decay faster than F^\ and 
that their decay exponents satisfy the inequality 

N(F {n) ) > iV(F (0) ) + nX n > 0. (96) 



The functions (|20[) that appear in the Gowdy metric (J3|) , expressed in 
terms of x = X\ and y = X2, are of the form (^), or in the case of the 
function A, which is defined by integral (|TTj) , the integrand is of the form 
(|93|). For the first three functions ( p0| ) the expansion is straightforward and 
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gives higher order corrections that are globally defined functions on the t6 
cylinder, whose decay exponents satisfy (|96|). 

The terms in the expansion of the function A are 

A< n > (t, 9) = 4 n) + a t n) o) dt> + f 4 n) (*> O^'> (97) 

./in J6»n 



00 

where a| n ^ and a.g are terms in the expansion of the functions at and a$ 

defined in (11) and (12), and A^ 1 a suitable constant. 

To determine the asymptotic behavior of the integrals note that, since 
both of 1 and ctjf* are ^-expandable functions we can integrate term by 
term. The time integral, when evaluated at the upper limit t, gives a £- 
expandable function with the decay exponent equal to the decay exponent 
of a[ n \ When evaluated at the lower limit to it gives a constant k(to,9o) 
that could spoil the decay properties of A^ n >. To prevent that we choose 



A 



(«) 



-k(t ,9 



Using (|9j|) the theta integral can be shown to be 



af\t,9')d9' 



< K {n) exp[-2 inf(l - c(0))nt], 



(98) 



for some constant > 0. 

We have chosen a special path of integration for A^ n \ but it is easy 
to show that A^ is path independent. Using the functional relation (|l4| ) 
between ag and at and reasoning similar to that that led to (|48|) we obtain 



a 



(«) 



(n) 



M - a t,e (99) 

which guarantees the path independence. To have each A^ defined globally 
on the t9 cylinder they must be periodic in 9, with period 2tt, which is 
equivalent to 



c (n) 



2~ 



a 



(n) 



(t,6)d9 = 0. 



(100) 



Equation ( |99| ) guarantees that C?W = 0, and we have to impose the con- 
straint at one time only. 

The decay exponent of is zero, which means that C(°) = will 
impose a constraint on the geodesic loop functions a, b, c and d. Evaluating 
C(°) and then going to the limit t — > oo we obtain 



2~ 



c(9) 



am 

b{9)) 



+ d'{9) 



d9 = 0. 



(101) 
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For n > 1 the decay exponents of are, according to (96), greater than 
1, which together with the constancy of (7'"' implies that = 0, without 
the need to impose any new constraints on a, b, c and d. 

So, we conclude that the expansion of the function A is globally de- 
fined provided the condition ( |101| ) is satisfied by the seed solution. Then, 
asymptotically for t — > oo, the zeroth order term in the expansion of A is 

A m (t , t _ i £ (102) 

and the higher order terms decay exponentially faster. 

Therefore, we have shown that for any geodesic loop spacetime with 
metric whose geodesic loop functions a, b > 0, c > and d satisfy the 



constraint ( 101 ) and 

c(0) < 1 or (103) 
c(0) arbitrary and a(9) + b(B) = const., (104) 

there exists a sequence <?( n ) of solutions of the perturbative Gowdy equations, 
to arbitrarily high order N, such that, asymptotically as t — > oo, 

AT 

£ 5 (n) - g {0) (105) 

n=0 



The functions a, b, c, and (i satisfying ( [L01D and ( |103| ) form an open sub- 
set in the set of all Gowdy initial data. This gives us a hope that the exact 
Gowdy solutions from an open subset Q of all Gowdy solutions asymptot- 
ically approach geodesic loop solutions. The geodesic loop approximation 
may not be valid for all Gowdy spacetimes however. Consider the speed 
of the exact Gowdy loops (pX|); there may exist Gowdy spacetimes with 
asymptotic speed lim^oo u = c > 1, and for such spacetimes our pertur- 
bative expansion is not asymptotically dominated by the zeroth order (the 
geodesic loop) term. In the case of c > 1, unless the condition (|104| ) is satis- 
fied, the higher order terms in the expansion exponentially increase instead 
of decrease compared to the zeroth order term, as t — » oo. 

All the rigorous results on Gowdy spacetimes obtained so far show that 
the geodesic loop approximation is indeed asymptotically valid. In the spe- 
cial case of the circular loop spacetimes, which because of circular symme- 
try can't satisfy the condition ( |l04j ), it was rigorously proven ||] that the 
geodesic loop approximation is asymptotically valid and that the asymp- 
totic speed c is always less than unity with all values between zero and 
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unity attained by some spacetime. Special condition ( |104j ) includes polar- 
ized Gowdy metrics which are rigorously known to asymptotically approach 
geodesic loop spacetimes with arbitrary large speed c [17]]. We will discuss 
these and some other exact results in the last section. 

For now, assume that our perturbative results are rigorously true, i.e. 
that the geodesic loop approximation is valid for all spacetimes with asymp- 
totic speed c < 1. The asymptotic behavior of curvature, as t — ► oo, in 
Gowdy spacetimes can be calculated using their corresponding geodesic loop 
asymptotes. The asymptotic behavior of the curvature in the geodesic loop 
spacetimes was analysed by Mansfield [11 1. He showed that for c / 1 the 



geodesic loop spacetime is curvature singular, as t — > oo, and 

R aMS R aPl S ~ \e- im (c{6? - l) 2 (c(9) 2 + 3)e< c W a+3 )*, (106) 

where 

m= f\ m m±mm M . (107) 

This means that the corresponding Gowdy spacetime with the asymptotic 
speed c < 1 is curvature singular, as t — > oo, and therefore inextendible 
beyond its maximal globally hyperbolic development, in agreement with 
the SCC. Mansfield also showed that for the special case of c = 1 and 
a + b =const. the spacetime is not curvature singular and is extendible. 



4 Concluding Remarks 

Chrusciel and Moncrief have derived some "light cone" and "higher-order- 
energy" estimates for Gowdy equations with a view towards rigorously prov- 
ing the geodesic loop asymptotic behavior suggested by the perturbation 
calculations presented here ||18|| . Several complications make such a proof 
significantly more difficult here than for polarized case treated in Refs.j|] 
and Ji7|]. Aside from the fact that now the basic equations are non-linear 
there is also the complication that not every geodesic loop can be realized 
as the asymptote of some exact solution. The reason for this is that our per- 
turbative calculations require c{6) < 1 to be applicable, unless some special 
non- generic condition such as a + b = const, is satisfied. 

Nevertheless, it is quite conceivable that exact field equations always 
force the asymptotic speed c below unity (unless the aforementioned special 
condition is satisfied) and thus force the solutions to achieve the geodesic 
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loop asymptotic behavior. An example of this phenomenon was found by 
Chrusciel and Moncrief who studied exact circular loop solutions, the motion 
of which is governed by a non-linear, second order ordinary differential equa- 
tion (for the "radius" of the loop) . The circularity is preserved by the exact 
field equations and such solutions are prevented, by their circular symme- 
try, from achieving the special condition which permits an asymptotic speed 
greater than or equal to unity (unless the geodesic loop behavior is violated). 
It was found, however, that every such circular loop solution does indeed 
asymptotically approach a geodesic loop (of speed less than unity) and that 
every value of asymptotic speed (strictly less than unity) is in fact achieved 
by some exact circular solution ||. 

The polarized solutions considered by Isenberg and Moncrief || satisfy 
the special condition mentioned above (a + b =const.) and can achieve arbi- 
trarily large asymptotic speeds while still realizing the asymptotic geodesic 
loop behavior. Not only is this asymptotic behavior universally satisfied 
by the polarized solutions but also one can rigorously justify computing 
the asymptotic behavior of the Riemann tensor or (for sufficiently smooth 
solutions) its covariant derivatives by means of computations made purely 
within the geodesic loop approximation. 

A further example of exact results which realize the geodesic loop asymp- 
totic behavior (this time with c = 1) has been found by Mansfield |ll]] who 
transformed a class of exact analytic solutions, derived using the Ernst for- 
malism, back to the Gowdy representation. This family of solutions (which 
is infinite dimensional and generically non-polarized) turns out to be none 
other than the set of "generalized Taub-NUT" spacetimes (restricted to the 
Gowdy symmetry class) defined on T 3 x R which develop compact Cauchy 
horizons instead of curvature singularity at the boundaries of their maxi- 
mal Cauchy development |14|, [15], . Much larger families of such general- 
ized Taub-NUT spacetimes (admitting generically only one spacelike Killing 
field) are defined in Refs. p3 and |l5j and have been used in Ref . J20|] as 
backgrounds for perturbative expansion. An infinite dimensional family of 
analytic, curvature singular solutions having only one Killing field (and de- 
fined on S 2 x S 1 x R) was constructed in by applying suitably chosen 
Geroch transformation to the generalized Taub-NUT solutions defined on 
S 3 x R. One hopes that these rigorously derived singular solutions will 
help one to understand the asymptotic behaviors of the perturbative so- 
lutions discussed in [ 20 1 . At our present level of understanding it seems 
not unreasonable to hope that the higher order perturbation methods dis- 
cussed here can be successfully applied to the study of completely general 
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(non-symmetric) cosmological solutions of Einstein's equations near their 
singular boundaries. 
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